Spin blocker made of semiconductor double quantum well using the Rashba effect 
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We propose a lateral spin-blockade device that uses an InGaAs/InAlAs double quantum well 
(DQW), where the values of the Rashba spin-orbit parameter cvr are opposite in sign but equal in 
magnitude between the constituent quantum wells (QW). By tuning the channel length of DQW 
and the magnitude of the externally applied in-plane magnetic field, one can block the transmission 
of one spin (e.g., spin-down) component, leading to a spin-polarized current. Such a spin-blocking 
effect, brought about by wave vector matching of the spin-split Fermi surfaces between the two QWs, 
paves the way for a new scheme of spin-polarized electric current generation for future spintronics 
applications based on semiconductor band engineering. 

PACS numbers: 72.25. Dc, 73.63.Hs, 85.75.-d 
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I. INTRODUCTION 

Semiconductor spintronics is a new paradigm for semi- 
conductor electronics which utilizes spins of electrons in 
addition to charges for device functionalities. Primary 
areas for potential applications include information stor- 
age, computing, and quantum information. One of the 
pedagogical devices in semiconductor spintronics is the 
spin field-effect transistor (spin-FET) proposed by Datta 
and DaSfi where a spin current injected from the source 
is controlled by gate via the Rashba spin-orbit interac- 
tion (S0I)f2,nl Three elemental technologies are required 
to realize switching in the spin-FET. The first is the 
generation and injection of spins into a semiconductor. 
The second is the manipulation of spins. The third is 
the detection of spins. Among these mechanisms, the 
generation of spin-polarized electrons without conven- 
tional ferromagnets has certain merits such as efficient 
spin injection and the absence of stray magnetic fields. 
While there have been many proposals to fulfill this 
purpose;^!^ the uses of the SOI provide the most popu- 
lar approachs. These include the intrinsic/extrinsic spin- 
Hall effect (SHE)ii"— spin-filtering devices using reson- 
nant tunnefing diodes (RTD)^"— those utilizing quan- 
tum point contact (QPC)^^"— and topologically pro- 
tected surface current in topological insulators (TI) J^iiS 
While the intrinsic mechanisms (intrinsic SHE and TI) 
potentially provide an ideal (the most efficient) source for 
the spin-polarized current, the actual utilization of these 
are hindered by the lack of their controllability. The ex- 
trinsic mechanisms, on the other hand, provide good con- 
trollability in the actual devices, whereas the efficiencies 
of the spin current generation are typically very low^S 

In this report, we propose a new spin-filtering device 
which provides both excellent controllability and effi- 
ciency, based on the extrinsic SOI mechanism, utilizing a 
narrow gap semiconductor double quantum well (DQW) 
struct ureJ^iSi"— In this device, the combination of the 
Rashba SOI and the external orbital magnetic field en- 



ables the perfect blockade of only one spin component, 
thereby generating a spin-polarized current. Our DQW- 
based device has advantages over the QPC-based devices 
in a sense that it is compatible with the conventional 
lithographic top-down processes and that a large (spin- 
polarized) current can be extracted from the device. 

We organize the paper as follows. We introduce the 
concept of the proposed device in the next section. This 
is followed by a description of tight-binding (TB) model 
that simulates the function of our device (Sec. Ill), where 
we also show how the parameter values for the TB model 
were chosen and that the resultant energy dispersion rela- 
tion indeed agreed with that of the effective mass model. 
Section IV is devoted to the discussions of the spin depen- 
dent conductance and the spin-filtering effect observed in 
the proposed device. Our conclusions are given in Sec. V. 



II. DEVICE CONCEPT 

A. Device description 

Figure l(a)(b) illustrates the proposed spin filter- 
ing device based on the DQW structure, where each 
QW (QWl and QW2), by itself, has an asymmetric 
confinement potential. The locally averaged electric 
fields within QWl and QW2, necessary for the Rashba 
SOI for each QW, are denoted as {Ez)i and {E^^, 
respectively!^ These electric fields satisfy the relation 
{Ez)i — — (£-2)2 > by symmetry where the z di- 
rection is normal to the sample surface. Such a DQW 
structure can be realized using (001) InP lattice-matched 
I1io.53Gao.47As/InQ.52Alo.48As material system as will be 
explained in Sec. IIIB. The thicknesses of QWl and QW2 
are both dqw ^ which are separated by a barrier layer with 
thickness de- The value of de should be so chosen that 
QWl and QW2 are only weakly coupled. The active part 
of the device, composed of DQW, has a length L in the x 
direction. Non-magnetic electrodes are attached to QWl 
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alone are invcrsionally asyminctric, which supports the 
presence of the Rashba splitting. However, if one sees 
the whole DQW as a single quantum mechanical system, 
it is inversionally symmetric, which is against the pres- 
ence of the Rashba splitting. Such dilemma is resolved 
as follows. If one sees QWl and QW2 as independent 
quantum wells, the Rashba spin splitting is indeed in- 
duced with a finite given wave vector k = (fcp, 0) in each 
QW. However, if one sees the DQW as a single quantum 
mechanical system, what we interpreted as the spin split- 
ting of a single quantum well above is now viewed as a 
spin degenerate subband splitting derived from a single 
DQW (see Figs. [2] and [S]). 



FIG. 1: (a) Schematic illustration of the proposed spin- 
blocking device made of double quantum well, (b) Sketch of 
the potential profile in the DQW region with zero magnetic 
field. The spin-dependent energy eigenvalues for the given 
wave vector value k = {kp, 0) due to the Rashba SOI are also 
indicated, (c) The equivalent tight-binding lattice model to 
simulate the transport properties of the electron along x. 




FIG. 2: (Golor online) Sketches of electron wave functions 
along z-direction at B = 0. (a) The case without the Rashba 
SOI or with k = (0, 0) for the in-plane wave vector, where we 
obtain the same orbital wave functions irrespective of spins, 
(b) The case with the Rashba SOI and k = (A;p,0). The 
Rashba splitting Ar is assumed to be much larger than the 
inter- well coupling fcoup. While the spin degeneracy is pre- 
served due to the inversion symmetry, the orbital wave func- 
tions are not identical between the different spins in each spin 
degenerate energy level. 



of the DQW device at both the left and right ends, be- 
tween which an electric current is passed through. The 
width of the device in the y direction is assumed to be 
much larger than L so that the periodic boundary con- 
dition is applicable. 

With the DQW structure described in Fig. [IJb), one 
may encounter the following dilemma. QWl and QW2 
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FIG. 3: (Golor online) (a) Illustration of unbiased DQW 
system without magnetic field, (b) A pair of spin degenerate 
eigenenergies for k = (fcp, 0) in the DQW. (c) Spin-dependent 
Fermi circles for independent QWl and QW2, where indicated 
by the arrows on the circles are the spin orientations, (d) 3D 
illustrations of the overlapped Fermi circles of the DQW with 
the given condition, where the vertical axis is z in the real 
space. 



B. Formation of the bonding and antibonding wave 

functions for a selected spin by the in-plane 

magnetic field at the Fermi circle points (±fcF,0) 



{H = Ho + Hr + Hz) reads 
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FIG. 4: (Color online) (a) Illustration of unbiased DQW in 
the presence of the magnetic field B = (0,5,0). (b) Spin- 
dependent eigenenergies in the DQW with k = (fcF,0) and 
B = Sac (see text for Bac). (c) Spin-dependent Fermi circles 
for independent QWl and QW2 at B = Sac, where we set z — 
at the middle of the center barrier layer, (d) 3D illustration 
of the overlapped Fermi circles of the DQW with the given 
condition. 



Let us consider the effect of applying an in-plane mag- 
netic field (B II y) to the DQW, which provides a key 
element to the proposed spin-blocking mechanism. Set- 
ting p — >■ p + eA in the Hamiltonian using the vector 
potential A = {Bz, 0,0), which leads to B = V x A = 
{0,B,0), the Hamiltonian of an electron in the DQW 



Hn = 



2mt: 



{k^ + eBz/hy + k^ 



and 



dz 2m* {z) dz 
Hr = a{z) [kyax - (k^ + eBz/h) ay] 
Hz = — 2.'7*/^B-Scrj/, 



V{z), 



(1) 



where Hq, H^. and Hz are the unperturbed Hamilto- 
nian, Rashba and Zeeman Hamiltonians, respectively, 
fiB = eh/2me is the Bohr magneton {me being the free 



'II 



electron mass), g* is the effective g- factor value, 
and m*{z) are the in-plane and out-of-plane effective 
masses, respectively, the latter incorporating different ef- 
fective mass values for the well and barrier materials. 
k = {kx, ky) is the in-plane wave vector of an electroni^ 
e is the elementary charge, h is Planck's constant di- 
vided by 27r. Ux and Uy are the Pauli spin matrices. 
For the Ino.53Gao.47As/Ino.52Alo.48As material system, 
we can safely ignore Hz and eBz/h relative to iJ^ and 
kx, respectively, and so we do hereafter j^ 

Using the material-dependent intrinsic constant ago 
for the Rashba SOI?^ a{z) can be written as 
(aso/e) {dV{z)/dz), where the contributions from the 
band discontinuities at the well-barrier interfaces are 
readily included in this formulation and should be ex- 
cluded in taking the derivative dV{z)/dz. We also de- 
fine the parameter value an, by ($qwi|q;(2:) |<i>Qwi) = 
aso(-Ez)i or equivalently - ($qw2| a(^;) |^qw2) = 
-aso {Ez)2i where |$qwi(QW2)) is the energy eigenstate 
of the independent QWl (QW2) in the confining direc- 
tion. 

In Eq. ([1]), we sec that applying a magnetic field 
(B II y) has the effect of shifting the Fermi circle in the 
kx direction by the magnitude —eBz/h. Setting the ori- 
gin of z at the middle of the barrier layer, the Fermi 
circles for QWl and QW2 are shifted oppositely along 
kx axis as shown in Fig. [D We then claim that the 
spin blocking effect is maximized when the magnitude 
of the Fermi circle shift becomes equal to the Rashba 
wave number ka, where ka = m,t,aR,/h'^. We call the cor- 
responding magnetic field as the "anticrossing" magnetic 
field i?ac = fika/e{z), where (z) = {^qw2\z\^qw2) = 
- ($Qwikl^QWi) > (see Fig.[ll). 

With this condition {B = Sac), defining the spin- up di- 
rection in y and neglecting the interaction between QWl 
and QW2, the front and back edges of the Fermi circles 
for spin-down electrons in the current direction (Hi) are 
at k = (fcF,0) and (— fcF,0), respectively, for both QWl 
and QW2. For spin-up electrons, those of QWl and QW2 
are at (i/cp + 2ka, 0) and (i/cp — 2ka, 0), respectively, as 
shown in Fig.[4l^c). If the interaction between QWl and 
QW2 is turned on, only the spin-down electrons form the 
bonding- and antibonding-like wave functions around the 
Fermi circle points (i/cp, 0) in the presence oi B = i?ac- 



Spin-blockade mechanism explained by the 
spin-selective flying qubit model 




FIG. 5: (Color online) (a) Energy dispersion relation of 
the proposed DQW when the in-plane magnetic field B = 
(0,-Bac,0) is applied, (b) Spin dependent trajectories of an 
electron which is injected in QWl from left, where the Fermi 
wave number matching condition is satisfied with the in-plane 
magnetic field B = (0, Bac, 0). 



Let us consider the Fermi circle points of the DQW on 
the kx axis in the presence of B = (0, i?ac,0). The k^ 
values of these points are — (kp + 2ka) (spin-up), — fcp 
(spin-down, doubly degenerate), — (fcp — S/cq.) (spin- up), 
fcp — 2ka (spin- up), fcp (spin-down, doubly degenerate), 
and /cp + 2ka (spin- up), if the inter- well coupling is ne- 
glected. The degeneracies at ±fcp are lifted if the inter- 
well coupling is turned on, i.e., fcp — >■ fcp ± fccoup with 
fccoup = tcoup{dE{kx) / dkx)~^ , where icoup and E{kx) are 
a half of the subband splitting energy and the energy 
dispersion relation of the unperturbed hamiltonian Hq 
with B = 0, respectively. We note that the states with 
±fcp — 2ka have a wave function along z nearly equal to 
|*QW2) (g) \ty), while that for ±fcp + 2ka is |$qwi) " 



Ity). The Fermi circle points (± (fcp -f- fcc 



, 0) and 



(± (fcp — /ccoup) , 0), on the other hand, form bonding- 
and antibonding-like wave functions between QWl and 
QW2, respectively. The superposition of these wave 
functions causes processional motion of electron between 
QWl and QW2 as we see below i^ 



$ 



Letting 
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|$QWi>® lly) =%l,0)®\iy) and 



QWl/ 

QW2) = l*QW2> ® \iy) = '(0,1) ® \iy) be the spin- 
down wave functions confined in QWl and QW2, re- 
spectively, the bonding- and antibonding-like wave func- 
tions, which are the spin-down energy eigenstates at 
(± (fcp -I- fccoup) , 0) and at (± (fcp — fccoup) 7 0) are, respec- 
tively. 
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Next, we consider how the wave function e*'^''^' 

in the left lead of the device propagates through the 
active part of the DQW device {0 < x < L). At a; = 0, 



$ 



U \ 
QWl/ - 72 



k{H') 



^i" 



from Eqs. ^ and ^, 



$ 



,-^y 



and 



$i" 



are 



where the plane wave parts for 
gi(feF+fecoup)a; j^jj(j gHfeF-fcooup)a:^ respectively, at the Fermi 
energy Ep. Thus, the wave function e^'^^^ *^qwi) ^-"^ 
the left lead {x < 0) connects to a superpositioned state 



*in.(-)>-^(^ 
1 COS yr^coupX ) 



mj 
,ik-px 



pi- recoup 3 



.ly \ 
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■i-y 



^ — ikc 



^i" 



V \\ _ 



$ 



QW2// 



""QWI / "*" ^^l'^ ("^coup2;) 

in the DQW (a; > 0) [the point indicated with 
in Fig. [5ja)]. This wave function portrays the pro- 
cessional motion of an spin-down electron between 
QWl and QW2 within the DQW, where the condi- 
tion L ~ Ln = {n ~ i) 7r/fccoup, n being an integer, 
makes an electron be backscattered at the end of 
the DQW within QW2. The backscattered wave 



function e" 



ik-p{x — L) 



$ 



QW2/ 



a,t X ~ L now con- 



nects to 



V2 



-, — ikc 



$^« ) _ e**^- 



^i" 



|-iSinffccoup2;'j ^QWl) +COS ffccoupS;'] *q'w2/| 

in the DQW [the point indicated with o in Fig. [SJa)], 
where x — x — L, which propagates back to QWl at 
a; = following the trajectory of "iff^- (x) ) backward as 
shown by the arrowed curves marked with in Fig.[5l^b). 
Electrons with up-spin, on the other hand, injected 
from the left lead to QWl, are transmitted straight to 



the right lead, i.e. 



ni, (^) 



— ^i{kp-\-2ka)x 



^QWi) as 



shown by the arrowed broken line in Fig. [SJb). 

The probabilities of finding spin-down and spin-up 
electrons in QWl at the position x = L upon inject- 
ing the same spin in QWl at position x ~ are given 
naively as 



T, 



T(QWi) 



$ 



QWl 



*a(i) 



(4) 



and 



T, 



;(QWi) 



$ 



QWl 



<"j(L))| = COS^ (fccoupi) (5) 



respectively. This means that the transmission probabil- 
ity of spin-down electron oscillates with the device length 
L, while that of spin-up electron is not (always unity). 
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FIG. 6: (Color online) (a) Equivalent circuit model for the 
gated DQW without the applied magnetic field, (b) Spin- 
dependent eigenenergies in the DQW with k = (fcF,0) and 
Vg = Ar/g [see (a) for Vg]. (c) Spin-dependent Fermi circle 
for independent QWl and QW2 with B = and Vg = Aa/e. 
(d) 3D illustration of the overlapped Fermi circles of the DQW 
with the given condition. 



D. Failure of the spin-orbit blockade by the Fermi 
circle matching with gate 

One may think that the same spin-blocking mechanism 
should be equally effective even if the size of the spin- 
dependent Fermi circles are controlled by the surface or 
back gate voltages, not by the in-plane magnetic field, 
but it is not the case as follows. Two Fermi circles of 
a selected spin can be made overlapped if the potential 
energies of QWl and QW2 arc shifted by a half of the 
Rashba splitting, i.e., — Ar/2 and Ar/2, respectively, or 
vice versa [sec Fig.[2][b) for ArJ. We note that the change 
in the value of or by this gating is almost negligible 
in the InGaAs/InAlAs DQW because the gate voltage 
required to shift the potential energies by ±Ar/2 is very 
small. It turned out, however, that the proposed spin- 
blockade mechanism can not be validated by gating at 
least within the scheme of the simple one-dimensional 
model presented here<^2i2i 
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FIG. 7: (Color online) (a) Energy dispersion relation of the 
gated DQW with 1/g = Ar/c [see Fig. [g a) for V^]. (b) Spin- 
dependent trajectories of an electrons which are injected in 
QWl from left, where the Fermi wave number matching con- 
dition is satisfied by gate {Vg = Ar/c). 



For example, with — Ap{/2 and Ar/2 for the shift of the 
potential energies in QWl and QW2, respectively, and 



with L = L„, the injected wave function e 



ikpx 



$: 



QWl/ 



at x = connects to 



^tn-i (^) 



which experiences a 
precessional motion between QWl and QW2 as in the 
case of the application of the in-plane magnetic field 
Bar [the point indicated with • in Fig. [7{a)]. How- 



ever, the wave function 
nects to e-»(feF-2fc„)(an-L: 



$^ 



\ 



QW2 / s-ftcr backscattcring in 
QW2 [O in Fig. Elja)], which propagates back to x' = 
through QW2, where it gets reflected again. Then, this 
wave function again experiences a precessional motion 
from QW2 to QWl and reach x = i in QWl, where 
the electron is ejected to the right lead [curves marked 
with in Fig. El^b)]. On the other hand, for a spin- up 



electron injected at x = 0, 



pi[k-p +2ka)x 



$ 



QWl 



simply 



propagates through QWl all the way to the right lead 



[broken line in Fig. [71[b)]. Therefore, both e'*^''^ 
and e'C'f +2fc, 



$ 



QWl/ 

*^QWi ) injected from the left lead of the 
device are eventually ejected to the right lead. This il- 
lustrates the failing mechanism of the spin-orbit blockade 
effect by gating instead of applying the in-plane magnetic 
field within the current simple one-dimensional modeli^ 



III. TIGHT-BINDING MODEL DESCRIPTION 
OF THE PROPOSED DEVICE 

A. Tight-binding model 

Keeping in mind that the proposed spin-filtering device 
has the translational symmetry along the y-dircction and 



the front edge of the Fermi circle in k^ carries the most 
of the electric current, we model the device by a tight- 
binding (TB) Hamiltonian^ consisting of a coupled one- 
dimensional chain as shown in Fig. 1(c), 



^ky - Yl 



^mky'^ma^ma' I / ^ ^ma'^jTam'^i^' <^' i 



^(To-' ;;t 



m,m' .(j,a"' 



where 



(6) 



be possible to design an efficient device whose width W is 
smaller than its length L to reduce the anti-spin-blocking 
effects while W is large enough to carry a large current. 
Once the Hamiltonian [Eq. ^] is set, the spin- 
dependent electric current in the right lead is given as 
Ii^^l ~ G^.xVsD in the linear response regime, where 
VsD is the spin-independent source-drain bias voltage. 
The zero-bias spin-dependent conductance G^,^ associ- 
ated with the right lead is expressed as follows in the 
Landauer-Biittiker formalism)^ 



vnkjj 



2mf, 



2to 



S^^' + 2atsoimz)ky [cra;]^^,(7) 



and 



X exp(=FJ0(m2)) 







^COUp'JfTfT' 



(m = m'zbe^j) 
(m = m' ± Bz) 
(otherwise) . 



Here, Cmo- {c\^^) is the annihilation (creation) operator 
of an electron at site m = {nix, mz) with spin a (=t,4,), 
where spin basis can be chosen in an arbitrary direction 
and rrix and ruz are the lattice site indices along the trans- 
port and out-of-plane directions, respectively. For m^ 
(=1 or 2), 1 and 2 represent QWl and QW2, respectively, 
to = h? /2rnT,a'^ is the orbital hopping parameter where a 
is the lattice spacing along the transport direction (a = 1 
nm in the present work). tso('7iz) = yn-z — §) aR/a is 
the Rashba SOI hopping parameter along the x-dircction 
within the mzth QW. The tunneling between QWl and 
QW2 is characterized by the inter-well coupling parame- 
ter tcoupj which is a half of the subband splitting for the 
unperturbed Hamiltonian Hq in the absence of B [see 
Eq. (1) and Fig.[5Ka)]. The effect of the applied magnetic 
field is incorporated as a form of Peierls phase factor 



(mz) 






(9) 



where $o = h/e and $ = 2Ba (z) are the magnetic flux 
quantum and the magnetic flux threading through the 
TB unit cell, respectively. The left and right leads are 
also described by the TB Hamiltonian Eq. © but the lat- 
tice site along the out-of-plane direction is restricted to 
"iz = 1 (QWl), and the SOI coupling and magnetic field 
are set to zero. While our TB Hamiltonian [Eq. (jH])] to- 
gether with Eq. pUj) below allows us to investigate on de- 
terioration of the spin coherence by finite ky component 
in the wave vector as following, our investigations in this 
article are restricted to the case with fcj, = to illustrate 
our fundamental idea of spin-blocking. First, the condi- 
tion of spin-blocking for ky = (e.g., B = Sac in Fig.[T 

is different from that for ky ^ Q and h 



± 



1,2 _ 7,2 
Zip /tj^. 



Second, the first term in H-r, in Eq. ([T]) [or the second 
term in Eq. ([7])] results in the precession of up- and down- 
spins in the y direction. The detailed evaluation of such 
anti-spin-blocking effects is being in progress. It should 
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ky (y'=\,i 



dET„ 



^AE) 



dfjE - Ef) 
dE 



(10) 

where Taa',ky {E) is a transmission probability of an elec- 
tron from the spin state a' {a' =t, i) in the left lead 
('g~j to the spin state a (cr =^"1, |) in the right lead, assuming 
the conservation of ky during transmission. In our simple 
model as mentioned above, we set ky = and skip the 
summation over ky. Then, T-^^.o(-E') = 2^it,o(^) ~ if we 
choose the spin basis in y, because Hk commutes with 
CTy. Thus, the spin- up and spin-down transmissions are 
treated separately. The spin-polarization of the electric 
current in the right lead is given by 



P = 



h - h 



L^ 



'i 



G^ — G 
G^ + G_ 



i 



i 



(11) 



The actual numerical values of T^a' ,k can be calculated 
using the Recursive Green Function technique from the 
TB Hamiltonian [Eq. dH)] as ellaboratcd in Ref. [H 



B. Values of the tight-binding parameters and the 
confirmation of the band dispersion of DQW 

The actual material system of our choice to make 
the proposed device is Ino.53Gao.47As/Ino.52Alo.48As 
DQW,-^ where the well width dgw for both QWl and 
QW2 is 10 nm and the barrier thickness d-Q ranges 1.5 
- 5.0 nm. The donor doping (Si^) in the carrier sup- 
plying layers above and below the DQW results in the 
electric fields along z axis within the DQW as explained 
in Sec. II. Other material specific parameters such as the 
energy band gap and the conduction band offset are given 
elsewhere.— 

The following values are used in our actual calcula- 
tions, to = h^ 12171*^^0^ = 0.81 eV from m*^^ = 0.047me 
and a ~ 1 nm. tgo ~ ctR/2a = 1.57 meV, where an = 

aso|e|WfV2 



2cseo 



3.14 X 10~^^ eVm assuming dV{z)/dz ~ 
at the center of the barrier layer between QWl and 
QW2 (z = 0). We note aso (intrinsic constant for 
the Rashba effect) = 25.28 A^ and es (dielectric con- 
stant) = 13.1 for Ino.53Gao.47As2^. The total sheet 
carrier density A^g°* was assumed to be 3.6x10^^ m~^ 
throughout the article. The values of tcoup (a half 
of the subband splitting) and (z) = ('I^qw2| -^ |'&qw2) 
are obtained from the self-consistent solutions of the 



one-band Poisson-Schrodinger equations of the unper- 
turbed Hamiltonian Hq with B = [Eq. ^] assum- 
ing the Neumann's boundary condition dV{z)/dz — 
at 2; = ± (^3/2 -I- 22) nm, which are just outside the 6 
nm thick, symmetrically placed carrier supplying layers 
above and below the DQW, where the donor density is 
1.8 X IQis m-2 (a half of 7V^°'). It turned out that the 
values of tcoup and (z) approximately obey the following 



phenomenological equations: tcou 
and (z)(nm) = -1.105 e-i^^' 
the unit for ds is nm. 



p(eV) = 11.66 6-109''=^ 
+ 5.744 + rfB/2, where 
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FIG. 8: (Color online) Energy eigenvalues in the symmetric 
DQW with de = 3 nm for a given wave vector k = (fcF,0) 
(fcp = 3.36 X 10* rn^^) as a function of the in-plane magnetic 
field, where the sohd and dashed curves correspond to the 
calculations with and without the inter- well coupling tcoup, 
respectively. Save stands for the averaged energy Save = {e\ + 
£i+£2+£2)/4- On top of the main figure, we depict the spin- 
dependent energy levels for tcoup = at three representative 
magnetic fields. 



The correspondence between the continuous effective 
mass model [Eq. (P)] and the TB model [Eq. ©] is 
made clear if one diagonalizes the following 4x4 TB 
Hamiltonian of the DQW part of the device in the 
spin(g)QW space for a unit cell, taking into account the 
/cj^-dependent Bloch phase assuming the translational 
symmetry along the x direction, which is equivalent to 
applying the periodic boundary condition. 
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are the couplings among the 4 sites in the spin(5 
space at the mxti^ lattice site and those between the rrixth 
and {nix ± l)th lattice sites, respectively. In Eq. (fT4)) . the 



matrix elements ig'^i a-^d ^0^2 ^^^ given by i^jn±e^ "^^ 
Eq. dll) with 771^ = 1 and 2 for QWl and QW2, respec- 
tively. 

Choosing the spin basis in y, iJspin(»QW j,^^-^ i-,^ block 
diagonalized and solved separately for each spin. The 
sorted spin-up and -down Hamiltonians are 



H^,l = 2to/Qw — 2-\/tQ + 1 



2 
so 



cos{kxa + (fi-f^i) 



''COUp 



tr 



cos{kxa- 0^,;) 



(15) 



in the QW space. 



= t. 



coup 



/2v*o+*so 



and 



\l 



± a 



We note icoup 

for t and — for |), where (f> ~ 
eBa{z)/h and a — tan^^ (tso/^o)- -^QW is the 2x2 
identity matrix in the QW space. The eigenvalues of this 
Hamiltonian, which provide the spin-dependent energy 
dispersion relations for the DQW system in the presence 
of the Rashba SOI and in-plane magnetic field, are 



ej'2 = 2io — 2-\/tQ -I- t^QCOskxacos(j)^^ 



i 



±\/^coup + 4 (^o + *io) ( sinfc,,a sinc?;^^ 



(16) 



where — and + signs are for subbands 1 and 2, respec- 
tively. One finds that these energies agree with the re- 
sults of the effective mass model in Eq. ([T]) within 1% 
at the Fermi energy by Taylor expanding the trigono- 
metric functions in Eq. (J16p . where the following values 
are sufficiently smaller than unity: k^a « k^a = 0.336, 
(/) « 4.3 X 10^3 (foj. 5 ^ 0.2 T and (z) = 7 nm) and 
a« tsoAo ^ 1-9 X 10-3. 

In Fig. in we plot the magnetic field dependence of 
these eigenenergies for de = 3 nm (<coup = 0.44 meV) and 
kx = kp = 3.36x 10^ m-i relative to their averaged value, 
together with the results when the coupling between 
QWl and QW2 is turned off (icoup = 0) . Here energy dif- 
ference between the spin degenerate pairs Ej ^nd e|' 



a.t B = 0, which is interpreted as the subband splitting in 
the symmetric DQW, is very close in value to the Rashba 
splitting energy 2aB.kF of independent QWl and QW2. 
This is because the orbital parts of the eigenfunctions for 
each spin-degenerate subband level are localized in either 
QWs depending on their spin state [Figl2][a)]. The slopes 
for B — )- T in Fig. [8] with icoup = are given exactly 
by zLhe (z) kp/mt. = ±5.97 meV/T. Upon increasing 

B, Fig. [5] shows that the energy levels ej and £3 anti- 
cross each other at i? = Sac = TJiT.aji / he (z) = 0.177 

T, where (j)i = (p — a = and Eq. ([T6| reduces to 
£ J 2 = 2io — 2-\/tQ -I- igfjcoskxa ± icoup- This means that 
the anticrossing energy gap found in Fig. [5] is exactly 
the subband splitting energy in the absence of both the 
Rashba SOI and the magnetic field [Fig. [2Ja)]. Thus the 
bonding and antibonding states are formed selectively 
for spin-down electrons. Spin-up wavefunctions are then 
localized in either QWl or QW2. As we further in- 
crease the magnetic field {B > i?ac), both the spin- up 
and -down wave functions are localized in either QWl or 
QW2 (sec the top pictures in Fig. [5]). 



IV. TRANSPORT ANALYSIS AND SPIN 
FILTERING PROPERTIES 
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In Fig. m we plotted the calculated values of the 
spin-dependent conductances G^^^ (t|| y) as a function 
of L at B == Sac for ds = 2 and 3 nm (T =5 K). 
We find that the change in the value of G^ as a func- 
tion of L is more pronounced than that of G^, where 
G^ even becomes zero at specific device lengths L*, n 
being an integer, while G^ varies only weakly with L. 
These qualitative behaviors are in agreement with the 
argument in Sec. II. We find that L\ = 648 nm and 
L*2 ~ 1946 nm for de = 2 nm, whereas L\ — 646 nm 
and L2 = 1939 nm in Sec. III. Similarly, L\ = 1914 
nm here agrees with L\ ~ 1906 nm in Sec. II for 
rfe = 3 nm. We note L„ = tt (?i — i) /fccoup, where 
^coup = t^oyip{dE{k) / dk)~^ « tcouprn*^\a/h^sin{kYa) in 
the TB model. 

Shown in Fig. [TU] are the magnetic field dependence of 
G^^l for ds = 2 nm and 3 nm, where the values of L 
are 646 and 1906 nm, respectively. For i? > 0, we ob- 
serve that G^ becomes zero &t B = 0.192 T and 0.177 T 
for de = 2 nm and 3 nm, respectively, i.e., at the anti- 
crossing magnetic field Sac- For i? < 0, on the other 
hand, G^ reaches zero at i? = — i?ac while G^ is kept 
close to unity, as is consistent with the symmetry found 
in the Hamiltonian [Eq. ([T])]. In Fig. [TOl we also rec- 
ognize that the variation of G^ (G^) around B = B^c 
(B = — -Bac) is more moderate for de = 2 nm than for 
de = 3 nm. This is because the bonding/antibonding 
states {1/V2{\^b) ± |*a)} ) at B = ±Bac are more ro- 
bust with the change of B for de = 2 nm than for de = 3 
nm because of the stronger inter- well coupling tcoup- 
The range of the magnetic field AB around B — i?ac 



FIG. 9: (Color online) Plots of the spin-dependent conduc- 
tances G-t-4, calculated for two different barrier thicknesses 

(a) de = 2 nm and (b) ds ~ 3 nm as a function of the device 
length L. The magnetic fields are fixed at their anticrossing 
values Bac = 0.192 T for (a) and 0.177 T for (b). The inset of 

(b) is the magnified view of the main panel between 260 nm 
and 300 nm for abscissa and between 0.9 e^ /h and 1.0 e^ /h for 
ordinate, where the period of the rapid oscillations is found 
to be Tr/fcp. 



within which the bonding/antibonding states persist can 
be roughly estimated by equating eAB (z) /h to kconp- 
We obtain AB = 0.24 T and 0.075 T for de = 2 nm and 
de = 3 nm, respectively, which agree well with the widths 
of conductance dips observed in Fig. 1101 Thus, measuring 
the width of the magneto-conductance dip would provide 
an estimate of the value of tcoup- 

We note that there are features in Fig. [9] that cannot 
be explained by the simple flying qubit model (Sec. II). 
The first is the rapid oscillation of the spin-dependent 
conductance as shown in the inset of Fig. Eljb). We no- 
tice that these oscillations (i) are more pronounced for 
shorter L, (ii) disappear as G^^^ approaches to zero or 
the maximum value e-//i. and (iii) damp away with in- 
creasing L. The period of the rapid oscillation, which is 
essentially constant with L, is found to be AL ~ 9.34 
nmjiSSi This is in good agreement with the value of ir/kp, 
which infers some resonance phenomena associated with 
the plane wave part of the electron wave function. De- 
tailed analysis using the multiple reflection model (Ap- 
pendix B) revealed that it is not the discontinuities in 
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FIG. 10: (Color online) Magnetic field dependence of the 
spin-dependent conductances for (a) ds ~ 2 nm and (b) da ~ 
3 nm. The device lengths for de = 2 nm and 3 nm are L = 646 
nm and 1906 nm, respectively. The anti-crossing magnetic 
fields Sac for ds = 2 nm and 3 nm are 0.192 T and 0.177 
T, respectively, as indicated in the figures. See text for the 
explanation of AB. 



the values of B and/or a^ in the Hamiltonian between 
the leads and the DQW but the finite probability am- 
plitude of the elctron wave function within QW2 which 
resulted from tunneling from QWl, that caused the wave 
function reflection at a; = and L. The reason for the 
disappearance of the rapid oscillations as G^ approaches 
to zero or e^/h is the disappearance of the wave function 
amplitude within QW2 at either x = or L with the 
condition L = L* or L = L* + 7r/2fccoup, respectively.— 
Second, the damping of the rapid oscillations with in- 
creasing L in Fig. [9] is due to the finite temperature as- 
sumed in our calculation. Considering the finite width 
Afc in the wave number values that participates in the 
electron transport, the criteria for having this rapid os- 
cillation is AkL < 1. Using Afc = kBT{dE/dky^ « 
8 X 10-*^ m"^ for T = 5 K, we obtain L < 1.25 nm, 
which explains the decay of the rapid oscillation with 
increasing L. 

Another feature in Fig. [9] that cannot be explained 
by the simple flying qubit model in Sec. II is the weak 
modulation of the spin-up conductance G^ with period 
7r/2fc^ in Fig. [SJ where k'^ is the corrected Rashba wave 



number (see below). This modulation of G^ can be 
explained qualitatively by the following extended qubit 
model (see Appendeix B for the exact model). It turned 
out that even the eigenfunctions of the spin-up electrons 
are more correctly described as superpositioned states 
between QWl and QW2 because the inter-well coupling 
icoup is not completely negligible for the spin-up elec- 
trons cither. These eigenfunctions still keep the features 
of antibonding-like or bonding-like wave functions in a 
sense whether they have a node or not, respectively. 



'i'D = 


= 5 


^QWl/ 


-Vl-5^ 


*QW2) . 
^QW2): 


(17) 


<) 


= Vi-s^ 


<J'qwi)+'^ 


(18) 



where < 6 < 1, apart from the plane wave part for the 
in-plane transport. We obtained 6 = 0.275 and 0.104 for 
de = 2 nm and 3 nm, respectively, from the eigenvec- 
tors of the spin- up Hamiltonian [Eq. (fT2)) ] . For electrons 
propagating in the positive x direction, the in-plane wave 
numbers for |$J^) and |$^) are kp + 2fc^ and kp — 2fc^, 
respectively, where the corrected Rashba wave number 
k'^ is ka + fc^oup/S^a incorporating the contribution of 
the inter-well coupling fcoup to the wave number shift 
(see Appendix A). After similar procedures as in Sec. II, 
which we call as the extended flying qubit model, we ob- 
tain the following expression for the transmission proba- 
bility of the spin-up electrons. 



T, 



T(QWi) 



1 



4(5^(1 



S^)sm^{2k'^L) 



(19) 



Thus, the periods of conductance modulation for spin-up 
electron are 7r/2fc^ = 677 nm and 792 nm for de = 2 nm 
and 3 nm, respectively, which are in good agreement with 
683 nm and 799 nm in the TB calculation (Fig.[9|). We 
note that the extended flying qubit model overestimates 
the oscillation amplitudes of G^, i.e., 4(5^(1-5^) = 0.2805 
and 0.0429 in units of e^/h for de = 2 nm and 3 nm, 
respectively, whereas the corresponding values in the TB 
calculations are about half of these. The discrepancy 
comes from the absence of both the multiple reflections 
of wave function between x = and L and the thermal 
averaging effect in our naive model. We derived the exact 
results assuming multiple reflections between x = and 
L which reproduced the TB results in Appendix B. 



V. CONCLUSION 

We proposed a lateral spin-blockade device using In- 
GaAs/InAlAs double quantum well (DQW), where the 
values of the Rashba spin-orbit parameter an are op- 
posite in sign but equal in magnitude between the con- 
stituent quantum wells (QW). The principle of the spin- 
blocking effect in the proposed device is in the spin- 
selective matching of the front and back edges of the 
spin-split Fermi circles [Fermi circle points at (±A:f,0)] 
between the two QWs, which is made possible by the 
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in-plane magnetic field B = (0, i?ac,0). The superposi- 
tion of the resulting bonding and antibonding wavefunc- 
tions that are formed for the selected (e.g., spin-down) 
electrons exclusively results in the precessional motion of 
electrons between the QWl and QW2, which is denoted 
as the "flying qubit" state. The "flying qubit" state can 
be blocked by depleting or etching away only the QW2 
part at a length of half-integer multiple of the precession 
wave length as depicted in Fig.ISjb). 

We would like to reiterate the features of the proposed 
spin-filtering device, (i) The proposed device is novel in 
a sense that it operates based on the conventional band 
theory and the Boltzmann transport theory, adding the 
spin degree of freedom in the form of the Rashba ef- 
fect. Thus, only the elementary level of quantum me- 
chanics and solid state physics is required to understand 
the principle of the device, (ii) In our simplified one- 
dimensional model where two singly-channeled leads are 
attached to the device, we obtain a perfect spin-blockade 
by the in-plane magnetic field B = (0, i?ac,0), whereas 
such spin-blocking fails to happen in the case of electrical 
control using the gate. The latter observation is found to 
be consistent with the recent theoretical results proven 
analytically for devices with two singly-channeled leads, 
(iii) The actual devices, however, will be prepared in the 
form of the two-dimensional electron gas (2DEG), where 
finite ky components (_L I) in the electron wave vector 
also participate in the electron transport. While some 
deterioration in the spin-polarization will be expected in 
the actual 2D model due to the finite ky components, 
the multi-channeled nature of the device would, in turn, 
allow the generation of spin-polarized current purely elec- 
trically. We can expect that the spin-polarized currents 
thus generated are fairly large due to the multi-channeled 
nature of the device, as compared to those generated by 
the QPC-based devices, for example, which open the pos- 
sibilities of future spintronics applications widely. 
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Appenxix A : Correction to the Rashba wave 
number ka due to the inter-w^ell coupling icoup 

The generalized Rashba wave number k'^ is associated 
with the energy difference ei^ ~ s\ with B = (0, Sac, 0) at 
the Fermi wave number as following [see Eq. (J16p ]. 



In the limit a — >■ 0, we have cf)^ ^ (f) + tan ^ (iso/*o) = 
2isoAo, where = eB^^al^z) jfi = tsoAo, iso/^o = 
/cqA, and tcoup — fi^kpkcoup/'n^u- Note to = fi^/2rnt,a?. 
Substituting these in Eq. ([20)) and letting a ^- 0, we have 



-t A 



Ah^k-pxka 



'II 



1,2 
coup 



4fc (i + ^'°"M ^^(^'-^ 



, (21) 



lij, — ftp 



for fccoup < ka- Comparing Eqs. (PH)) and (PTj) . we obtain 



ka ~ ka 



coup 
Ska ■ 



(22) 



Appendix B: Multiple reflection model 

Consier a DQW device with length L. Let t„im' 
(jn, m! = 1, 2) be the nonnnal quantum mechanical trans- 
mission amplitude of electron from x = in QW?n' to 
X = L in QWtti. Similarly we let rmm' be the nomi- 
nal transmission amplitude of electron from x = L in 
QWr?T,' to a; = in QWrn. The multiple reflection model 
states that the overall trasmission amplitude from QWl 
at a; = to QWl at a; = L is given by 



it. 



til + il2'"22i21 + il2'"22i22?'22i21 + ' ' 
in + il2'"22 |l + (i227'22) + (i22''22)' 
il2?'22i21 



iii + 



1 - i22?'22 ' 



•}t21 

(23) 



where the transmission probability Ttot is given by jttot | ■ 
The values of tmm' and rmm' per spin are obtained by 
the generalized flying qubit model as following. 

Choosing the spin quantization axis in y, the eigen- 
functions in DQW without spin are e.±<>'v+^^)=^ |<j)^^) and 
^±i(k^-Ak)x 1^^^ at E = Ey and B = (0,5ac,0), where 



|$b) = Vl-<^2|$QWi)+'5|*QW2), (24) 

|$a) = <5|$QWl)-v/T^^|$QW2) (25) 



(0 < (5 < 1). We note Afc = fccoup or 2fc^ for spin-up or 
spin-down electrons, respectively, and '(Vl — <^^, 5) and 
*((5, — Vl — 5^) arc the eigenvectors of Ec^. ([T5|) . We solve 
these for |$qwi) and |$qw2)- 



e\ = 



2\/i?oup + 4(to+*so) (sinfc-pasin^^ 



^ 4fc' ^^(^-) 



dkrr 



(20) 



I^QWl) = VT^^|$b)+'5|$a) (26) 

|*QW2) = <5|$b>-v/l-'52|$a> (27) 

We let these wave functions propagete from a; = to 



11 



QW2 



QWl 







FIG. 11: Schematic illustration for symmetry consideration. 
The nominal transmission amplitude through QW2 with neg- 
ative fca; (r22 in the left panel) is found to be equal to that 
though QWl with positive kx (in in the right panel). 



X = L. 



|$QWi> (atx^O) 



Ak^L 



|g,AfcL^^ - <52 1$^) + e-'^'^'^S |$a)} 



Ak-pL \ ( AAkL 



e-'- [{e'^''^ - 2iS^aiTi{AkL)} |$qwi; 



+2iS^/l - (52sin {AkL) |$qw2) 
(at X — L) 



|*QW2) (ata;==0) 



AkpL 



{e'^^^5 |$b) + e-*^^^ ^Jl-5-^ |$a> } 



Ak-pL 



2i5y/l~Psin{l\kL) |$qwi) 



+ {' 



-iAkL 



2i6hm{AkL)}\<^QW2) 



(at X ^ L) 



(28) 



From these we obtain 



til = e'^''^ {e'^''^ ~2i5''-siVL{AkL)] , (29) 
ii2 = 6*'=''^ {2i(5v/l - <52 sin(AfcL)} , (30) 

e'*^"-^ {2i5 Vl - 5^ sin(AfcL)| (31) 



t21 

and 



t22 = e'''''^ {e-"^''^ + 2iS^sm{AkL)} . (32) 

Now we let |$qw2) propagate from x = L to a; = to 
obtain r22- After a simple symmetry consideration in 
Fig.im we notice that it is equal to in. 

r22 = hi^e''"'^{e'^''^~2iS^sm{AkL)}. (33) 



Substituting Eqs. 
have 



and ^ into Eq. 



AkpL , 



g.,.,^ ^^ _ g2ifcpL^ ^giAfeL _ -^2 siii(fcp7^)] 
1 - e2*feFirii 



, (34) 



where Tn = |fii|^ = 1 - 4(5^(1 - (5^) sin2(A/i:L) is the 
result of the extended qubit model. We obtain the overall 
transmission probability as 



l^totP 



2 [1 - cos(2fcFL)] Tn 



(35) 



1 - 2rii cos(2fcFL) + Ti^i ' 
We confirmed that this reproduces the TB results quan- 
titatively. We note that the spin-dependence of Ttot is 
realized only through Tn. Thermal average of Ttot is 
also straightforward. 
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